A Semi-Lagrangian Finite Volume Method for Newtonian Contraction Flows

A new finite volume method for solving the incompressible Navier–Stokes equations is presented. The main features of this method are the location of the velocity components and pressure on different staggered grids and a semi-Lagrangian method for the treatment of convection. An interpolation procedure based on area-weighting is used for the convection part of the computation. The method is applied to flow through a constricted channel, and results are obtained for Reynolds numbers, based on half the flow rate, up to 1000. The behavior of the vortex in the salient corner is investigated qualitatively and quantitatively, and excellent agreement is found with the numerical results of Dennis and Smith [Proc. Roy. Soc. London A, 372 (1980), pp. 393–414] and the asymptotic theory of Smith [J. Fluid Mech., 90 (1979), pp. 725–754].